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Abstract 

The effective action and the mean fermion number density of graphene in 
constant external magnetic field at finite temperature and density are calculated. 
Closed expressions for these are given and their variation with temperature are 
studied. It is found that the mean fermion number density peaks around a par- 
ticular temperature, depending on the chemical potential at low temperatures. 
This feature is interpreted as 'condensation of fermions ipip^ in graphene. In 
future, it is interesting to extend and explore this calculation and the work of 
the reference [20] for the case of Graphyne [23]. 



Graphene is a honeycomb lattice of Carbon atoms with remarkable elec- 
tronic properties [1,2]. Many of these properties are a low energy phenomena, 
governed by massless Dirac equation for the electrons in the tight-binding (hop- 
ping) model. This has been discovered in a seminal contribution of Semenoff [3] 
and exemplified by Mecklenberg and Regan [4] and Ando [5]. Thus, the exci- 
tations of graphene in the tight-binding approximation near the Fermi surface 
are described by (2-|-l)-dimensional massless Dirac equation , with velocities 
vp ^ 10^ ms~^. This description is accurate theoretically [3], [6], neglecting 
many body effects and experimentally [7,8], as the quasiparticles in graphene 
exhibit the linear dispersion relation E = pvp (as if they are massless relativistic 
particles of speed vp) and by the observation of a relativistic analogue of inte- 
ger quantum Hall effect. Recently, it is shown mathematically in [20], that the 
linear dispersion relation (Dirac points) is generic in (non-relativistic) quantum 
crystals having the Honeycomb lattice symmetry. It is important to realize that 
the quasi-particles of graphene have to be described by two-component wave 
functions, necessary to define relative contributions of sublattices [9]. Such a 
two-component description is similar to spinor wave functions in QED, with 
the 'spin index' for graphene indicating sublattices rather than the real spin of 
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electrons. 



A graphenc layer has been extracted out of graphite [10] and such a discov- 
ery of 2-dimensional system attracted tremendous interest in theoretical physics, 
such as massless Dirac fermion, quantum Hall effect and even fractional quan- 
tum Hall effect [11], providing a framework to use the methods developed in 
quantum gaiige field theories in graphene. Many of graphene physics is de- 
scribed by the Quantum Field Theory (QFT) [21]. It is remarkable to note that 
many High Energy Physics (HEP) Experiments can be realised with the graphene 
table-top experiments. This statements is written in the spirit of Moore 's Law 
for Electronics. It is proposed in [22] how certain experimental results for carbon 
nanotubes (very closely related to graphene) could be re-interpreted in terms of 
the Kaluza-Klein mode of extra dimension. Now the physics is more interesting 
with the advent of the new material Graphyne [23j, which is more versatile than 
graphene due to its direction- dependent dirac cones. A functional determinant 
approach with zeta function regularization has been employed by Beneventano 
and Santangelo [12] to study the quantum Hall effect in graphene and by Ben- 
eventano, Giacconi, Santangelo and Soldatic [13] to study the one-loop effects 
for massless Dirac fields coupled to electromagnetic backgrounds both at zero 
and finite temperature and density. The present author have studied SU{2) 
Yang-Mills theory at finite temperature in (3-|-l)-dimensions for understanding 
the vacuum in Quantum Chromo Dynamics [14] and so it is worthwhile to use 
the methods developed in [14] in graphene. 

It is the purpose of this contribution to derive closed expression for the 
effective action Fe// for graphene (quasiparticles in graphene) in an external 
constant magnetic field at finite temperature and density. The 'mean fermionic 
number density' Np = ^ ^^f^ i where /U is the chemical potential, is calculated. 
We have used the e-regularization method of Salam and Strathdee [15] to reg- 
ularize the divergences to obtain a finite action. The behaviour of the effective 
action and Np with temperature has not been reported so far in the literature 
to the best of our knowledge. 

Our results for T^ff and Np are plotted against temperature for represen- 
tative values of /i. We find that the effective action peaks around a particular 
temperature for a given (jl. The 'mean fermionic number density' also peaks 
around the same temperature. This peaking of Np is interpreted as a possible 
'fermion pair condensation' or 'quasi-particle pair condensation' in graphene. 



We consider massless electron in graphene coupled to an external magnetic 
field {U{1) gauge field) in (2-|-l)-dimensions, described by the Lagrangian den- 
sity 

C = ijijt'D^^; A={-t^,0,Bx), (1) 
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where wc use h = vp = 1 units and the eonstant magnetic field is perpendicular 
to the plane of graphene. /x is the chemical potential for the electrons. For the 
Euclidean version of (1) (useful to study the theory at finite temperature) we 
use 

7° = -(^3 ; 7^ = ; 7^ = ; = 26'^, 

p = iP-e^. (2) 

The eigenvalue equation 

{i^-e4)tp = XTp, (3) 

gives 

A = ±v/2A^eB + (w„ + M)2, (4) 

where the harmonic oscillator quantum number A'' = 0,l,2,---oo and the Mat- 
subara frequency u;,, = (2n + 1)^; n = — ooto +oo, with (3 = k being the 
Boltzmann constant. Equation (4) gives the Landau levels of the fermions in 
graphene. Prom the (Euclidean) partition function 



j [#][#] 



-J d?x Ci 



the effective action Fe// = tog Z = log det{ip — e^) is given by 

^^/^ = "4^(y S ^ iog(j^{2NeB+{u;r.+f^n), (5) 

''^ JV=0 n=-oo 

where the prefactor (^ff^ is the harmonic oscillator degeneracy factor and a 
dimensionful parameter A is introduced to render the argument of the logarithm 
dimensionlciss. In odd dimensions (such as here), there is another inequivalent 
representation for the gamma matrices, obtained for example, by reversing the 
sign of one or all the three gamma matrices in (2). Inclusion of this introduces 
a factor two for Fe//. 

Now we use the e-regularization of Salam and Strathdee [15] (see also [16]) 
to write 



A''=0 n=— oo 



in the limit e — )• 0. After a Wick rotation t — >■ —it and carrying out the A/'-sum, 
we find 



4/3\2ttJ ^ \er(e)J Jo (1 - e-2^«■B/A2^ 

g-^(a;„+Mf_ (7) 
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The n-sum can be done using Jacobi ^3-function to give 

4l3\2TrJ\er{e)J Jo (1 - e-2"B/A2) 



Using the property of ^3-function [17], namely 



e3{z,it) = ^3(|,^), (9) 



we find 



8^\2TTJ\eT{e)J Jo (l _ g-s^^-B/A^) 
Rewriting Os{z^t) as 

00 

^3(2, t) = ^ e*""' e^""^ (11) 
£=-00 

equation (10) becomes 

'^^^ ~ 8V^\27r J \eTie)J Jo (1 - g-^-s/A^) 

00 ,222 

{1 + 2^ e~^^ cos(£7r + £/i/3)}. (12) 
^=1 



In (12), the term without the £-sum has the integral over r is evaluated [18] 
and the limit e — >■ is taken to give 

where the superscript denotes the contribution to Te// coming from the term 
without the ^-sum. This contribution corresponds to the effective action at zero 

temperature and zero chemical potential, and agrees with [13] and [19]. In (12), 
the term involving the f-sum is evaluated, first by writing 

00 00 

E 2neB ^— v 2neB 
e = 1 + e 

n=0 1 



1 

(1 _e-2eBr/A2) 
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and then using the function K^{xz) [18] and taking the limit e ^ 0, to give 

^S, = i(f)i:-(^' + WH|r(l) 

+ t '-^^^^^ K.HPV^)}, (14) 



n=l 

SO that the effective becomes (adding (13) and (14)) 



1 feB\ ^ cos(£7r + £/x/3) ^ 

+ ^[^) 2. «7 



+ ^ 2^2-^/ l3iV2neBKi {pN2neB)). (15) 

n=l 

Finally, using Ki{z) = [18], the effective action (15) is written as 

+ ' f; '""' + f:^--"-^}). (16) 

€=1 n=l 

This is one of our main results, giving the effective energy density of the quasi- 
particles in graphene in a constant magnetic field at finite temperature and 
chemical potential. It is important to notice that the dimensionful parameter 
A, introduced in (5), has disappeared, unlike in (3+l)-dimensional gauge field 
theories [14,16]. The 'mean fermion number density' can be found using Np = 

1 OTVi 

In order to give the variation of Fg// with temperature, it is convenient to 
introduce 



X = /3^/2iB; y = /x/? = xy, (17) 



so that 



and 

Np 1 'A sin(^7r + Ixy) 
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= -TT^I-:^ tan 

^ g g sin(^7r + fay) g-^y^r^i^ ^^g^ 

^=1 n=l ^ 



In Figure 1, — — ^ is plotted against i = j^eB '^'^^ V — 0-5 {[x 

/2eB 



(*)■ 



0.5v2ei?). The effective energy density peaks around ^ = 0.18. In Figure 
2, —, — is shown against for y = 0.5. The 'mean fermion number 

density' peaks around — 0.18. In Figures 3 and 4, we give them for y — 0.8. 
We have found for y = 0.4, 0.5, 0.8, 1.0, the peaking occurs at different values of 
x~^ . The minimum value of y for the peak to occur is found to be y = 0.32. 



This feature of peaking of Np is interpreted as 'condensation of quasi-particle 
charge density'. Since the fermion number or charge density in this framework 
is proportional to it is the pair ■0?/; (like Cooper pair) that condenses at low 
temperature and hence no violation of Pauli principle. This feature is different 
from the chiral symmetry breaking as it occurs around a particular temperature. 
Second, chiral symmetry breaking occurs when ^^^|m=o 7^ Oj where m is the 
fermion mass. Since the role of /i is similar to that of m, we see from (18) and 
(19) that I^^o = and so we realize a genuine condensation phenomenon. 

For the further deeper interpretation of the condensation we can look [24] and 
try to link our results in terms of material sciences. 
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Figure 1: Variation of scaled effective energy density with scaled temperature 
for y=0.5 
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Figure 2: Variation of scaled charge density with scaled temperature for y=0.5 
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Figure 3: Variation of scaled effective energy density with scaled temperature 
for y=0.8 
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Figure 4: Variation of scaled charge density with scaled temperature for y=0.8 
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